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Buoyant magnetic flux ropes in a magnetized stellar envelope 

Idealized numerical 2.5-D MHD simulations 
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ABSTRACT 



Context. The context of this paper is buoyant toroidal magnetic flux ropes, which is a part of flux tube dynamo theory and the 
framework of solar-like magnetic activity. 

Aims. The aim is to investigate how twisted magnetic flux ropes interact with a simple magnetized stellar model envelope — a 
magnetic "convection zone" — especially to examine how the twisted magnetic field component of a flux rope interacts with a 
poloidal magnetic field in the convection zone. 

Methods. Both the flux ropes and the atmosphere are modelled as idealized 2. 5- dimensional concepts using high resolution 
numerical magneto-hydrodynamic (MHD) simulations. 

Results. It is illustrated that twisted toroidal magnetic flux ropes can interact with a poloidal magnetic field in the atmosphere 
to cause a change in both the buoyant rise dynamics and the flux rope's geometrical shape. The details of these changes depend 
primarily on the polarity and strength of the atmospheric field relative to the field strength of the flux rope. It is suggested 
that the effects could be verified observationally. 
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1. Introduction 

Buoyant magnetic flux tubes are an essential part of 
the framework of current theories of dynamo action in 
both the Sun and solar-like stars: it is widely believed 
that flux tubes are formed by some combination of ro- 
tational shear and turbulent convection near the bot- 
tom of the convection zones (CZ) of these stars. When 
sufficiently buoyant the magnetic flux rise in the form 
of tubular (toroidal) il-shaped loops of magnetic field 
lines. Rising under the infiuence of rotational forces, 
the loops finally emerge as slightly asymmetric and 
tilted bipolar magnetic regions at the surface, see e.g. 
Fan et al. (1994)| and [Caligari et al. (1995)| During the 
last decade, it has been established that if they do ex- 
ist, these fiux tubes must in fact be flux ropes of in- 
tertwined or twisted field lines; otherwise they would 
quickly be disrupted by a magnetic "mushrooming" in- 



ficiently high, corresponding to a critical value of the field 
line pitch angle '^c approximately determined by equality 
between the energy density of the twisted field component 
and the ram pressure due to the rise of the flux rope e.g. 



Moreno-Insertis & Emonet (1996) 



stability, see e.g. Emonet & Moreno-Insertis (1998) and 



A fundamental ingredient in all theories of solar-like 
dynamos is a process where a relatively weak poloidal 
magnetic field is turned into a toroidal field by differen- 
tial rotation: it is this toroidal fleld that eventually gives 
rise to the flux ropes, which at the end of their existence 
return a poloidal fleld component to the CZ, thereby com- 
pleting the magnetic activity cycle. This dynamo process 
is gradual and different cycles overlap. 

Buoyant magnetic flux ropes have been exten- 
sively study through numerical simulations, cf. the re- 
Dorch (2002), and [Fan (2004)[ and the most 
high resolution study in two dimensions by 



views m 
recent 



Cheung et al. (2006) However, so far studies of the inter- 



Dorch fc Nordlund ( 1998)| It has been shown that the in- 
stability is inhibited if the degree of systematic twist is suf- 
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action of flux ropes with poloidal fields have solely dealt 
with flux ropes that emerge into the solar corona, e.g. 
Archontis et al. (2004)| i.e. not with a magnetic CZ. 
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The question addressed in this paper is how the po- 
larity and strength of the fleld in a magnetized convec- 
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tion zone may affect the rise and evolution of twisted 
flux ropes. When the predominantly toroidal buoyant flux 
ropes rise through the CZ, they encounter a poloidal mag- 
netic field that has a component perpendicular to the 
ropes' axes. However, since the fiux ropes are twisted, 
the transversal field components may be either parallel 
or anti-parallel to the average magnetic field in the CZ, 

1. e. the field in the two-dimensional plane perpendicular 
their toroidal axis. 

Several questions relevant to the theory of flux tube 
dynamos then emerge: e.g. is the rise faster or slower and 
do more or less magnetic fiux reach the surface in the 
presence of a magnetized convection zone; how does a 
magnetized convection zone affect the required amount 
of twist. In this paper I attempt to shed light on these 
questions by presenting results from idealized numerical 
2.5-dimensional MHD simulations of the cross-sections of 
buoyant fiux ropes that interact with a horizontal mag- 
netic layer while they rise. 

2. Model 



The full compressible MHD-equations are solved 
ing the stagger-code by Galsgaard and others. 



Galsgaard fc Nordlund (1997) 



dp 

dt 

dt 
de 
dt 
dB 
'dt 



= -V-pu, 

= -VP+j X B - V • (puu) 

= -V-(eu)-FV-u + gv-f 

= V X (u X B) +7/V^B, 



us- 
cf. 



(1) 
(2) 
(3) 
(4) 



Here p is the fluid density, u is the velocity, P the gas 
pressure, j is the electric current density, B is the magnetic 
field density and e is the internal energy. In Eq. ||2J) j x B 
is the Lorentz force and Qv and Qj are the viscous and 
Joule dissipation respectively. 

The equations are solved numerically on a staggered 
mesh using derivatives and interpolations that are of 6th 
and 5th order respectively in a numerical scheme that con- 
serves V-B = exactly. The time stepping is implemented 
by a third order predictor-corrector method. Numerical so- 
lutions are obtained on a high-resolution two-dimensional 
Cartesian grid of 2048^ points. The initial set-ups of the 
models are twofold, consisting of a snapshot of a stratified 
and magnetized adiabatic CZ-model, and of an idealized 
twisted magnetic fiux rope. 

Initially the entropy in the interior of the ropes are set 
equal to that in the external medium leading to a buoy- 
ancy equal to 1/7/? (with 7 = 5/3). The coordinate system 
is chosen such that x is the vertical coordinate, y the hor- 
izontal coordinate and z is the coordinate along the axis 
which is parallel to the rope. 

The fiux ropes move in a polytropic atmosphere P = 
Po{pI PgV J where Pq and po are the quantities at the initial 
position of the flux rope corresponding to a pressure scale 



height of Hpo = 59 Mm. The top pressure scale height 
is one tenth of the bottom scale height HpQ. The upper 
boundary is well below the position of an imagined photo- 
sphere. Horizontally the boundaries are periodic and the 
physical size of the computational box is 136 Mm in both 
the vertical and horizontal dimensions corresponding to a 
height of 2.3 Hpg. The initial position of the flux ropes are 
set to a; = 0.27 HpQ. 

The initial twist of the flux ropes are given by 



B^ = a{r/Rf B,e-('-/«)' 



(5) 
(6) 



where B^ is the parallel (axial) and B^ the transversal 
component of the magnetic held with respect to the ropes' 
main axes. Bp is the amplitude of the fleld, R the ra- 
dius, and a is a field line pitch parameter. The critical 
pitch of a twisted fiux rope needed to prevent the dis- 
rupting instability depends on the ratio of the rope's ra- 
dius to the local pressure scale height: for thin flux ropes 
this ratio is small and in the case of the thin fiux ropes 
studied by Emonet fc Moreno-Insertis (1998)[ the critical 
pitch angle was 10°. In the following I set the initial ra- 
dius to i?o = 0.14 Hpo, i.e. the tubes are non-thin yielding 
a critical angle of $c ~ 25° according to the expression 



by Moreno-Insertis fc Emonet (1996) To be on the safe 
side, a = 1.75 is chosen, so that the initial maximum field 
line pitch angle is 35.6°, occurring at the border of the 
rope at r = ■\/3/2 R for the topology given in Eq. ©. 
Incidentally this is low enough that a three-dimensional 
fiux rope would not be kink unstable, since to achieve a 
substantial growth rate this requires a pitch exceeding 45°, 



e.g. Galsgaard fc Nordlund (1997) 



Table 1 . List of numerical simulations and their basic param- 
eters: model name, ratio of CZ field to the rope's twisted field 
component e, sign of the latter ratio, and ratio of the rope's 
axial field to the CZ field strength x (only the value of Bcz 
vary). The horizontal dividing line indicates the approximate 
position of the lower solar limit, i.e. x© ~ 10. 



Model 


e 


II 


X 








N/A 


N/A 


lA 


0.044 




90 


IB 


0.044 


+ 


90 


2A 


0.22 




18 


2B 


0.22 


+ 


18 


3A 


0.44 




9 


3B 


0.44 


+ 


9 


4A 


2.2 




1.8 


4B 


2.2 


+ 


1.8 



With a grid of 2048 points in each dimension, the fiux 
ropes are resolved by more than 200 grid points while 
the resolution of the pressure scale height Hpo is closer 
to 1000 grid points. From a computational point of view 
one has to set the ropes' plasma /3's lower than the solar 
values (10^ — 10'') to reduce the computational time scale 
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to a reasonable value: Pq = 107 is chosen for the initial 
strength of the buoyancy on the ropes' axes. 

The magnetic layer is implemented by adding a mag- 
netic field of By — Bqz to the hydrostatic equilibrium at 
heights X > 0.9 Hpo above the bottom of the computa- 
tional box. Numerically this results in a steep gradient at 
that height (see the discussion on diffusion in Section 4). 

For convenience I furthermore define x = Bj-ope/ Bcz 
for the ratio of the ropes' axial fields to the CZ field 
strength, and e = Bqz/B^ for the ratio of CZ field to a 
characteristic value of the ropes' twisted field components, 
e.g. its maximum value. In the simulations presented here 
both e and x &re varied by a factor of 50 between the 
weakest and strongest CZ field strength. For all models 
ex has the same value w 4. An estimated value of x for 
the Sun is given by the ratio of the canonical toroidal field 
strength of 100 kC, to a poloidal field of at most equipar- 
tition Boq ~ 10 kC, i.e. xo > 10- 

It is trivial to calculate a simple limit for the strength 
of the envelope field Bcz relative to the strength of rope 
by assuming a force balance between the buoyancy of the 
flux rope and the oppositely directed tension caused by a 
dent of magnitude R in field lines in the CZ field layer: 
with an initial buoyancy of 1/7/3 the order of magnitude 
of X should exceed 



2 Hpo 



7- 



R 



(7) 



if buoyancy is to overcome the magnetic tension from the 
field lines in the layer. The size of the ropes in this model 
results in xo ~ 5. For the Sun, the approximate lower 
bound fulfills xq > Xc and hence solar flux ropes are ex- 
pected to be able to rise through the poloidal field layer 
as assumed by the flux tube dynamo. 

3. Results 

Several models were run, varying parameters such as nu- 
merical resolution and Bcz, but not the ropes' initial field 
strengths. Tabled lists the models that are discussed in 
this paper. Except from models 4A and 4B all models have 
X > Xc- 

The nine models consist of eight models of flux ropes 
rising into a magnetized layer corresponding to a poloidal 
field in a CZ and a single reference model where no CZ 
field was present model 0. 

3.1. Flux conservation and numerical dissipation 

With a weak or absent envelope magnetic field, the total 
magnetic flux through the plane perpendicular to a rope's 
axis 4>z remains constant: it only decreases by 0.1 promille 
during the simulation runs and hence there is virtually no 
flux loss due to (unwanted) numerical dissipation. The 
total magnetic flux through the horizontal boundaries is 
another matter: the total flux $y through the horizontal 
boundary decreases constantly by typically 0.5 promille 
per time unit, leading to a loss of 1-3% during an entire 



run. This loss could in theory be due to magnetic field 
being pushed through the upper boundary. However, the 
upper boundary remains closed to the passage of plasma 
and as will be discussed subsequently, the flux loss results 
primarily from natural physical diffusion of the CZ mag- 
netic field layer. 

3.2. Buoyant rise of the ropes 

The single simulation (model 0, lacking a magnetized 
stratification) is completely identical to that presented 
in the review by Dorch (2002) and elsewhere (albeit here 
in much higher numerical resolution): as the rope rises 
and expands it approaches a "terminal rise phase" where 
it rises with a speed approaching a terminal velocity 
Vt (determined by the force balance between buoyancy 
and drag). In this phase the rope e.g. oscillates due to 
its differential buoyancy and rise more or less in the 
same way as an adiabatically rising, non-stretching flux 
tube obeying an internal polytropic equation of state, see 



Dorch fc Nordlund (19981 and|Dorch (^002) 



For all the models in Tab.^ the initial behavior of the 
ropes is unaffected by the presence of the overlying mag- 
netized layer. When a rope begins to rise due to its buoy- 
ancy however, it eventually approaches the layer: while 
the plasma in front of the rope moves away to "let the 
rope through" , part of the mass — in front of the rope (be- 
cause of mass continuity, Eq. ^ — moves upwards along 
with the rope and pushes on the horizontal magnetic field 
lines, thereby causing an indent of compressed field in the 
layer. At some point the tension in the dented field lines 
becomes sufficient to withstand the compression by the 
upward fiow, and the forefront of the rope and the nearly 
horizontal field lines begin to approach each other faster. 
In the models, this happens at a time of t ~ 50 Tg, with 

The most pronounced effect is clearly the case where 
the CZ field is the strongest, corresponding to e = 2.2 
in the simulations model 4A and model 4B that have an 
envelope field of the same order of magnitude as the ini- 
tial strength of the ffux rope and hence X < Xc- The 
rise of these ropes is quickly halted and they never get 
past the threshold to the magnetized envelope: they re- 
main below the layer while performing damped oscil- 
lations eventually becoming a perpendicular (toroidal) 
magnetic layer, similarly to what happens to weak fiux 
tubes in photospherical simulations, e.g 
and |Archontis et al. (2004)[ 



Magara (2001) 



The simulation sets numbered 3A-I-3B, 2A-I-2B and 
lA-l-lB respectively, are more interesting in the solar con- 
text. These simulations correspond to situations where the 
rope's axial strength is approximately 10, 20 and 90 times 
that in the magnetized envelope: this covers an order of 
magnitude in the ratio of the envelope's field to the rope's 
twisted field component, e. Figure ^ compares B^ for the 
two first of these simulation sets to the reference model 
at a late time {t ~ 95 Tg): the effects of the presence of 
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Fig. 1. Grey-scale panel showing snapshots of Bz at t ^ 32.6 time units for models with moderately magnetized 
envelopes. Minimum (zero) and maximum magnetic field strengths are stated for every snapshot (in units of 354 kG) . 
Vertical dimension corresponds to x in the simulations. Shown are model (shown twice), and in the upper row models 
2 A and 2B (initial Bqz/B^ = 18) and in the lower row models 3A and 3B (initial Bcz/B^ = 9). Also shown are 
horizontal lines indicating the initial position of the magnetized layer. 




2G' 40 60 80 100 120 140 20 40 60 80 100 120 140 

Time (t/rj Time (t/rj 



Fig. 2. Left: the height of the flux ropes (in terms of Hpo) as a function of time for six models indicated by symbols. 
In order of increasing strength of the magnetic layer the models are: lA (squares), IB (crosses), 2A (stars), 2B (plus), 
3 A (diamonds) and 3B (triangles). Also shown is reference model (solid line). The vertical dashed line indicates the 
initial lower boundary of the magnetized layer in the CZ. Right: speed of the flux ropes as a function of time (in units 
of the sound speed at the initial position of the ropes). Here the thick solid line illustrates the theoretically expected 
asymptotic behavior of the rise speed for an un-magnetized CZ. 



the poloidal layer are clear, and these will be discussed in Figure El (left) shows the height of the flux ropes as a 
the following. function of time: in the case of e = 0.22 and e — 0.44 the 
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CZ field lowers the speed of rise, compared to the reference 
model. Neither of the two models with e = 0.44 reach 
very far before being halted within the magnetized layer, 
most easily seen in case of model 3B in Fig.|21that has an 
envelope-parallel field, where the flux rope is trapped at a 
height of ~ 1.1 Hpo, just in side the magnetic layer. 

Figure 12 (right) shows the same trend in terms of the 
rise speed: the expected evolution for a flux tube rising in 
a field-free environment is an asymptotical approach to a 
terminal rise speed. However, with a magnetic CZ corre- 
sponding to the high e values, after reaching a maximum, 
the rise speed plummets to zero, when the ropes encounter 
the magnetized layer at time t ~ 50 Tg. An exception may 
be model 2A that has an anti-parallel twist field: at the 
end of the simulation this rope still had a substantial non- 
zero speed, although declining. 

On the one hand, as one may naively expect in the four 
situations just discussed, the parallel-field layers lower the 
speed more than anti-parallel layers does. This is because 
of the impossibility of reconnection between the two mag- 
netic flux systems in this case, that would otherwise lower 
field line tension by removing transversal fleld from the 
region above the apex of the flux ropes. On the other 
hand, consider the cases with e = 0.044 and x = 90: these 
ropes actually rise marginally faster than the CZ field- 
free reference model, while again the simulation with an 
anti-parallel twist rises the fastest of the two (see Fig.E)). 

Figure El (left) shows the decrease of the axial mag- 
netic field strength relative to the field of the reference 
model as a function of time. The axial field of the reference 
model decreases in a manner similarly to an adiabatically 
expanding one-dimensional tube, see e.g. Dorch (2002) 
For the models with non-zero Bqz, the deviation from 
the reference model is substantial: for higher e the field 
decreases more slowly with time than the simple model. 
This is because in the high e models the flux ropes rise 
more slowly and hence expand more slowly. In situations 
where the flux ropes are halted, it is no longer the time 
scale of the rise that dominates the decreasing field, but 
the diffusion time, which is much longer. In case of low e 
the rise is slightly faster than the field-free reference case 
and then so are the ropes' expansion and the decrease of 
the field strength implying negative — -Bref as seen in 

Fig.m 

In addition, the ropes that have an anti-parallel 
twisted field component have lower axial field strengths 
than their parallel-field counterparts when they begin to 
rise. The explanation of this lies in what happens to the 
protective twist that surrounds the flux ropes' cores (see 
the discussion below). 

4. Reconnecting the twist 

Figure 121 (right) shows the flux through the xz-plane (i.e. 
through the plane spanned by the flux rope's axis and the 
vertical) as a function of time: at first one may find it 
surprising that it seems that the low e models loose flux 
more rapidly than the ones with a stronger magnetized 



layer. However, there are two effects that may cause a de- 
clining y-flux: flux loss out of the computational domain 
through advection and magnetic dissipation (diffusion). 
The magnetic layer could in theory be expelled from the 
computational domain via advection, but the average ver- 
tical transport velocity on the top boundary remains neg- 
ligible, because of the closed upper boundary condition. 

The second agent for reducing the y-flux is diffusion, 
but reconnection between the twisted fleld of the flux rope 
and the magnetic layer only changes the field lines' con- 
nectivity and does not contribute to diffusion of the y-flux. 
In stead, diffusion takes place at the bottom of the mag- 
netic layer where the vertical gradient of the field initially 
is very sharp. Hence, natural physical diffusion causes the 
initial behavior of to be identical for all the models up 
to around time equals 40-50 r^: the flux decay is propor- 
tional to slow exponential decay on a time scale of roughly 

6000 Tg. 

After the ropes begin to interact with the magnetic 
layer, the results again depend on e: there is a notice- 
able difference between the high and low e models. As 
mentioned, in the latter models begin to decay more 
rapidly, as the magnetic layer is mixed into the lower ly- 
ing atmosphere, which initially was field-free: this mix- 
ing introduces further diffusion of the y-flux through new 
magnetic gradients developing at small scales. 

Figure21(left) shows the evolution of the fleld line pitch 
angle ^max for the two models with e = 0.044. These two 
models have oppositely directed twisted fleld components 
compared to the direction of the envelope field. The pitch 
clearly develops in a manner different from the field-free 
reference model, and the evolution depends on the direc- 
tion of the twist: as expected the pitch of the rope in model 
IB increases beyond the reference model, as the parallel 
field lines of the rope's twisted field is aligned up against 
the envelope field, which has the same polarity just in front 
of the rope. In the opposite case the pitch decreases when 
the rope in model lA encounters the magnetized layer: the 
core protecting twist is gradually peeled way by field line 
reconnection, revealing the less pitched field lines belong- 
ing to the 'interior' of the rope. It is this effect that causes 
the dependence on the sign of the twisted rope component 
on e.g. the dechning y-fiux ^y (Fig. Q and the decrease 
of the axial field strength B^ (Fig.jJl right). 

FigureQl (right) shows a section around the apex of the 
flux rope of model 2A: the magnetic field vectors change 
direction across the current sheet that surrounds the up- 
per part of the flux rope. In Fig. [3 (left), the field lines 
corresponding to the vectors in Fig. ^ are shown for a 
larger section near the rope's apex. On this scale the be- 
havior of the magnetic field is very reminiscent of recon- 
necting field lines around a magnetic null point, leading 
to an X-type topology: near the central vertical line above 
the rope's center, the upper almost stationary field meets 
an incoming oppositely directed field, leading to recon- 
nection across the null point which separates the two field 
line systems. The reconnected highly bend field lines move 
away almost horizontally. However, in Fig. (right), the 
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60 80 100 120 140 20 40 60 80 100 120 

Time (t/rj Time (t/r,) 



Fig. 3. Left: maximum axial magnetic field strength as function of time (not height). The field strength for the same 
six models as in Fig. [5] are plotted with the field strength of the reference model subtracted. Right: the magnetic flux 
through the xz-plane as function of time for the same six models as in Fig. [51 The vertical dashed line indicates 
approximately the time when the ropes' centers pass the initial lower boundary of the magnetized layer. 




Time (t/rj Distance from rope center (Mm) 



Fig. 4. Left: the maximum field line pitch angles are shown as a function of time for the models lA (squares) and IB 
(crosses) that correspond to e = 0.044 and oppositely directed twisted field components. Also shown is the reference 
model (stars). Right: subsection of the plane around the apex of the rope in model 2A. Shown arc the magnetic field 
(vectors) and current density (grey scale). 



view is zoomed in on the current sheet by a factor of 10: 
the reconnecting region of the field clearly does not have 
a simple X-type topology, but rather a more complicated 
situation arises where closed magnetic field lines (plas- 
moids) are formed in the center of the current sheet, e.g. 



Shibata et al. (1995) and Archontis et al. (2006) 



The interaction of the flux rope with the overlying 
magnetic layer not only changes the flux ropes' twist 
through field line reconnection, but also alters the geo- 
metrical shape of the ropes: the general tendency is that 
ropes with a twist directed parallel to the CZ field become 
more roundish as they rise into and become encapsulated 
in the magnetic layer, see e.g. the model 2B in Fig. ^ — 
ropes with anti-parallel twist loose some of their protective 
transversal 'shielding' and their axial fields are much more 
steeply inclined at the ropes' upper boundaries (e.g. model 
2A in Fig.P). 



5. Discussion and conclusion 

It can be argued that two-dimensional models have several 
restrictions compared to the full three-dimensional case: 
e.g. there is a tendency for enhanced fiows due to the lack 
of the extra degrees of freedom associated with three di- 
mensions. Therefore, one can consider the questions posed 
in this paper to be addressed here only to a certain order. 

On the one hand, the initial naive assumption that one 
could have — that the effect of the magnetized layer would 
be to destroy a flux rope with a twist anti-parallel to the 
polarity of the layer — is wrong. In fact, the opposite is 
true: even though reconnection within the magnetic layer 
reduces the twist of the flux ropes, it is easier for them to 
penetrate the poloidal layer, because the reconnection of 
field lines allows the ropes to "carve" their way through 
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Fig. 5. Magnetic field lines near the front of the flux rope in model 2A at time — 20. The upper most field lines belong 
to the magnetic layer in the CZ and point to the right, whereas the lower more curved field lines from the flux rope 
points towards the left. Left: view of the rope's apex on a scale of ~ 25 Mm. Right: zoom at the apex where the field 
lines are nearly horizontal. 



the layer. This is provided, however, that the flux ropes 
contain enough transversal flux. 

On the other hand, the simple estimate of Xc ~ 5 in Eq. 
|7|seems to hold in the simulations: in the models with x < 
Xc, the ropes' buoyancy are strongly damped. However, 
while ropes with x > Xc a-re allowed to rise, they may still 
be halted in the CZ, but the more correct critical limit 
will have to be a function also of twist, i.e. Xc = xi^)- 

It has been shown here that the flux ropes' dynamics 
depend on both the strength and sign of the poloidal field: 
for a relatively strong poloidal field the rise is slower and 
therefore the flux ropes reach a lower height in the same 
amount of time. In certain cases the rise can be completely 
halted. It turns out that anti-parallel twisted ropes reach 
higher and faster. When it comes to the flux ropes' topol- 
ogy, the geometrical shape of ropes also depends on the 
strength and sign of the poloidal field: the apex of anti- 
parallel twisted ropes are flatter and have steeper mag- 
netic gradients in their axial field. Ideally this could be 
observed as a seemingly faster passage through horizontal 
layers (emergence) of the axial part of the fiux rope, rela- 
tive to cases with more moderate vertical gradients in the 



axial field: e.g. Archontis et al. (2004) have shown that in 
fact a steep axial gradient is required by flux emergence 
(their Eq. 10). As to the flux ropes' twisted field com- 
ponents, reconnection reduces the twist of anti-parallel 
ropes. If ropes on either sides of the equator are oppo- 
sitely twisted, this could then in principle be falsified ob- 
servationally by comparing observations of poloidal field 
strength in emerging northern and southern active regions. 

The most fundamental problems remaining are those 
of the origin of the twist, and the question of how it arises. 



One may speculate that twisted field lines could be gener- 
ated in large-scale flux bundles located near the bottom of 
the convection zone connecting across the solar equator: 
such flux bundles would experience a rotating motion since 
their lower parts are located in a region rotating slower 
than their uppermost parts. This rotation would twist the 
magnetic field lines and this twist could be transmitted to 
the parts of the flux bundle at slightly higher latitudes, 
thereby possibly giving rise to a twisted toroidal flux sys- 
tem. The sign of a twisted component generated in this 
geometric manner would always be pointing in the same 
direction (northwards in the Sun). Hence, in consecutive 
11-year half-cycles of the solar dynamo, the twist would 
alternate between being anti-parallel and parallel to the 
poloidal field (from which the corresponding toroidal field 
was generated in the respective cycle). Such an alterna- 
tion would result a 22-year cycle in the amount of twist in 
the emerging parts of the toroidal flux ropes, e.g. in the 
helicity of bipolar magnetic regions. 

Fully three-dimensional MHD simulations extending 
the prese nt study are under wa y, using the numerical 
scheme of 



Archontis et al. (2004) 
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